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Abstract 

We discuss on the critical exponent of rj/s for a fluid, and propose a new exponent-less 
measure of fluidity based on a mode-mode coupling theory. This exhibits a remarkable 
universality for fluids obeying a liquid-gas phase transition both in hydrodynamic as 
well as in nonhydrodynamic region. We show that this result is independent of the 
choice of the fluid dynamics, viz., relativistic or nonrelativistic. Quark-Gluon Plasma, 
being a hot relativistic and a nearly perfect fluid produced in relativistic heavy-ion 
collisions, is expected to obey the same universality constrained by both the viscous 
and the thermal flow modes in it. We also show that if the elliptic flow data in RHIC 
puts a constraint on rj /s then the new fluidity measure for Quark-Gluon Plasma in turn 
also restricts the other transport coefficient, viz., the thermal conductivity. 
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Recent results from the Relativistic Heavy Ion Collider (RHIC) at BNL reveal 
surprising and intriguing dynamical properties of the Quark-Gluon Plasma (QGP). In 
non-central collisions, the anisotropy with respect to the reaction plane, i.e., the elliptic 
flow coefficient V21|21j can well be described up to transverse momenta of order 1.5 
GeV/c by the nearly ideal hydrodynamics with small shear viscosity, rj . This is much 
smaller than that obtained in the perturbative quark-gluon plasma liSJ. This suggests 
that the matter produced in the early phase of the RHIC collisions is strongly inter- 
acting with nearly perfect fluidity represented by a very small ratio of shear viscosity 
to entropy density, rj/s. The supporting a[;guments for it have come from various di- 
rections: viz, viscous hydrodynamics 1^]^ sets an upper bound, rj/s < {5/4-7Z); the 
gauge-gravity dual theory [6] conjectures at universal lower bound, rj/s > I /An; and 
the Heisenberg's uncertainty principle implies a lower bound L^. JJ, rj/s ^ \ (we use 
units with Ti = ks = c = \). The fluidity of QCD fluid has also been computed from 
numerical simulation in lattice |8], non-perturbative models |0,[l3l and has been com- 
pared with commonly known fluids. Moreover, for fluids rj is either a 
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nondiverging or a very weakly diverging quantity ifTil ^ with an exponent —0.04. 
Thus, the fluidity of a hot viscous QCD matter remains an open as well as interesting 
problem that requires a much desired exploration of the fluid mechanics. 

One can now ask what is the physical quantity that would reflect a good measure of 
fluidity of a hot viscous fluid or in this context how much relevant is the rj /s, which is 
formed by two quantities, completely different in nature as Tj is a dynamical transport 
coefficient whereas s is a thermodynamic one. A possible explanation lies within the 
general structure of any diffusion constant may it be thermal diffusion, Si^ related to 
the diffusive decay of longitudinal components of momentum fluctuation or the viscous 
diffusion, Sl-q related to the diffusive decay of transverse components of momentum 
fluctuation of fluids. In this letter we try to address this important aspect of a hot fluid 
based on fluid dynamics. 

Dynamical properties of a many particle system can be investigated by employing 
an external probe, which disturbs the system only slightly in its equilibrium state, and 
by measuring the response of the system to this external perturbation. A large number 
of experiments belong to this category, such as studies of various line shapes, acoustic 
attenuation, and transport behaviour. In all these experiments, one probes the dynam- 
ical behaviour of the spontaneous fluctuations in the equilibrium state. In general, the 
fluctuations are related to the correlation function, which provide important inputs for 
quantitative calculations of complicated many-body system. 

The density-density correlation function is defined as 



(T„„(r,f) - (5n(r,f)5n(0,0)) , (1) 

where the angular bracket denotes an equilibrium ensemble average, and 5n(r,f) = 
n(r,f) — (n), is the local deviation of the number density n(r,f) from the equilibrium 
value of the number density. The spectral function can be obtained by Fourier transfor- 
mation of ([T]i as 

P«„(q, (O) - [d\ r dt (J„„(r,f) , (2) 



and the static correlation function can be defined as 

Tnn(q) j ^p„„(q,«) . (3) 

Now, within hydrodynamical description the correlation function can be calculated 
only when wave number of the dynamical density fluctuation, q is much smaller that 
the inverse correlation length, ^ (q^ << 1), or, equivalently wavelength, A of density 
fluctuation is appreciably larger than the correlation length (A >> ^). In the case of 
nonrelativistic fluids, the dynamical density fluctuation based on Navier-Stokes equa- 
tion f]3\ in the hydrodynamic limit as well in the nonhydrodynamic limit (viz., around 



critical point) has been studied Ill6ll for Newtonian fluids in details. Recently, a rel- 
ativistic generalisation of density fluctuation is obtained [17] using dissipative rela- 
tivistic fluid dynamics coupled with conserved quantities and the dynamical structure 



2 



function is 



P«n(q,«)/T„„(q) = 



7/ ©2 + (^^^2)2 y 



(4) 



where 7 is the ratio of specific heats, Cp — T{dS/dT)p at constant pressure to C„ = 
T{dS/dT)„ at constant density with S is the equihbrium entropy per particle, is 
the speed of sound propagation with the damping constant and is the thermal 
diffusivity. The higher-order terms involving the quantities ^k-^/c? and S^sq/cs are 
neglectedQ. The thermal diffusivity is defined as, 

(5) 



noCp 

where K is the thermal conductivity and «o is the equilibrium number density of fluid. 
Note that the nonrelativistic equivalence of noCp in (|5]) is poCp, where po is the equi- 
librium mass density. 

The sound wave damping constant is given as 

&s^^Aj-l) + ^^^+c^j'-(—-2S)^ap), (6) 

Wo \Wo ) 

where Tj and 1^ are, respectively, shear and bulk viscosities, wq is the equilibrium en- 
thalpy density and ap is the thermal expansivity at constant pressure. The third term in 
^ is purely relativistic correction whereas the second term is the minimal relativistic 
one due to the appearance of wq instead of po. 

The dynamical structure function of density fluctuation for relativistic dissipative 
fluids in (|4]i is a sum of three Lorentzian line shapes with general form /(o)) = 2r/ [F^ + 
(o — (b')2], centered at the frequency ta' with a half-width at half maximum given by 
r. The Rayleigh component, is centered about = 0, with a half width 

r« = ^^<?' . (7) 

The Brillouin doublets, are located symmetrically about = at the frequencies O)^ = 
ztcjQ', each one with a half-width given by 

'^B^\%q^ ■ (8) 



A relativistic first order equation, sucli as that by Landau or by Eckart, is parabolic and formally violates 
the causality, and hence acausal. The causality problem can be circumvented with the inclusion of the second- 
order in derivative expansion fl flTlfTill . However, this does not affect the spectral function of dynamical 
density fluctuation. 
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Figure 1 : (Colour online) The structure function p„„ / T„„ as a function of (O in the minimal relativistic fluid 
(red curve), nonrelativistic fluid (black curve) and their difference (blue curve) for parameters fTJ q = 
0.1 fm"', chemical potential /J = 200 MeV, T = 200 MeV, rj/(Ho5) = t;/{"oS) = 0.3 and KT/{noS) = 0.6. 
The minimal relativistic correction shows up in the Brillouin doublets but not in the Rayleigh component 
centered at O) = (see text). 

We note that in ^ the transverse component of the momentum fluctuation decou- 
ples from the density fluctuation but admits a diffusive solution that relaxes exponen- 
tially with the viscous diffusion, !^ri as 



which reduces to rj/po in the nonrelativistic limit. The knowledge of the line widths 
in ^ and (O along with the three diffusion equations (|5]), ^ and ^ are sufficient to 
determine the three transport coefficients, K, Tj and ^. 

Now we note that the peak and the width of the Brillouin doublets obviously change 
due to the relativistic corrections of the sound mode, f^, in In contrast the Rayleigh 
peak as well as the width Fr remain same for both nonrelativistic and relativistic fluids 
because the thermal diffusion does not change due to the equivalence between 
noCp and poCp. We demonstrate this aspect in Fig. [T] with only the minimal relativistic 
corrections in (|6]l. Further, the line width Fr is controlled by the dominant behaviour of 
either k or Cp in The behaviour of will control that of Fb through the relative 
behaviours of k/C„, rj and ^. Also the strength of the various peaks in (01) depends 
on y. This hydrodynamics predictions have been used in the limit T to study the 
behaviour of relativistic fluid (v7\ around critical point. The fundamental assumption 
of hydrodynamic, << 1, close to Tc is no longer valid but it was argued that for a 
given q there will be always a temperature range over which hydrodynamics predictions 
should be reliable. Around the critical point 7 diverges faster than the viscosities (viz., 
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one can easily find from (|4|l that the sound mode gets attenuated and the only 
surviving soft mode is the diffusive thermal mode. This has important consequence as 
we will see below. 

Generally, the long wavelength part of the density fluctuation around critical point 
is very intense that could induce a velocity gradient in the boundary of the fluid, which 
would dissipate energy. This dissipation can be interpreted in terms of the mode-mode 
coupling fl9[|20ll . For instance, as discussed above, the heat flow mode (thermal diffu- 
sion) decays into shear viscous mode plus a thermal mode. One can obtain this decay 
mode in terms of the relevant transport coefficients. 

Information on critical behaviour of fluids are generally extracted from scatter- 
ing experiments performed on the fluids [2M- To analyse these data Swinney and 



Henry 112211 . based on mode-mode-coupling theory, obtained a particular dimensionless 
combination of measured quantities for wide range of fluids obeying liquid-gas phase 
transition as 

r^^. (10) 

Using the measured values of F/;, Tj and the correlation length, £, for different temper- 
atures and scattering angles (q — 4nsm{0/2)/X, 9 is the angle of scattering), obtained 
for various sample of multiple component fluids [21], when plotted F* as a function 
of q£, all fluids obeying liquid-gas phase transition fall on a single universal curve as 
shown in Fig. |2] The single curve describes the critical behaviour not only along the 
critical isochore and the coexistence curve, but also along any thermodynamic path in 
the critical region. It is clear from Fig. |2]that, 

/ ^ for<7(^<l (hydrodynamic domain), ^^^^ 
\ 1 for > 1 (near critical point), 

provided ^^j^f^ ~ 1 , and using ^ it becomes 

TJ 1 

'^'"^^(i?.^)-'r = 6^' ^^^^ 

which is independent of the nature of the fluid despite the various microscopic details. 
This is also consistent with Stokes formula for the mobility of a sphere of radius ^ 
moving through a viscous liquid with a nondiverging shear viscosity, rj . 

Now, the equality, F* = 1 in (fTTT i basically holds at and near the critical point, 
indicating the length scale matching around the critical point. This implies that instead 
of diverging without any upper bound, the correlation length of fluid becomes of the 
order of the wavelength of the density fluctuation (^ ^ A). On the other hand this 
also suggests that in together is a nondiverging quantity around critical point 



since both of them are expected to have same critical exponent 11151] but opposite in 
sign. The dimensionless quantity, ,^riK associated with the viscous and the thermal 
modes together, could then be regarded as a good fluidity measure of a hot viscous 
fluid obeying liquid-gas phase transition, which remains at a universal value l/(67r) 
for both nonhydrodynamical as well as in hydrodynamical region. 
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Figure 2: F* as a function of for several several substances that cover a wide range of molar mass, 
chemical structure and complexity and it is adopted from (22]. Note that the measured data are from various 
experiments and the detailed references of those can also be found in 

We now note that the important quantity that enters in ( fTTl i and (fT2b is the Rayleigh 
width, Yr. The analysis, based on the Navier-Stokes theory presented in the beginning 
of this letter, clearly exhibits that the relativistic correction does not show up in the 
Rayleigh width Yr in It indicates that the fluidity, ,j^r\K is, obviously, independent 
of the choice of the (viz., relativistic or nonrelativistic) fluid dynamics. This possibly 
indicates that a good fluid, either relativistic or nonrelativistic, is always a good fluid 
irrespective of the other details of the microscopic degrees of freedom. 

In a heavy-ion collisions, the produced matter has a strong flow initially along the 
beam axis. If the matter is viscous the viscosity then counteracts this by reducing the 
effective longitudinal pressure and thus enhancing it in the transverse plane. Therefore, 
the viscosity leads to a larger radial flow than that of the ideal fluid and also to an 
angular modulated radial flow (V2) due to the velocity gradient. This is indeed the 
situation in the experiments at RHIC thatprobably probed a QGP at a temperature 
little above Tq with a small value of T] /i ||3, S [13, ElD implying that the QCD matter 
produced in such experiments is a nearly perfect fluid and likely to obey a liquid-gas 
phase transition with a critical end point The spectral analysis for relativistic 

fluid in (|4|i or in Fig. [T]represents that a small density perturbation propagates at the 
speed of sound and is eventually damped near the critical point. The only surviving 
mode is the diffusive heat flow mode. The observed fluidity could be interpreted on the 
basis of mode-mode coupling theory as the decay of a longitudinal heat flow mode in 
QGP to a viscous and a thermal mode in which the transverse viscous mode reduces 
the thermal flow in the QGP. The QGP fluid produced in RHIC belongs to the same 
universality class defined by the fluidity measure in (fT2] l that interrelates the relevant 
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transport coefficients, viz., the shear viscosity rj, and the thermal conductivity K. 

We now recall that for fluids 77 is either be nondiverging or weakly diverging (with 
an exponent ~ —0.04) |7, 14]. One can expect that 77 /s should also be a nondiverging 
or weakly diverging quantity as entropy generally does not show any critical exponent. 
Surprisingly, for simple fluids as well as for QGP rj/s (as obtained from lattice QCD 
simulation |8] and also from the viscous hydrodynamic calculations ||3>El[l2l in view 
of recent RHIC data [2]) shows an exponent (see Fig. 3 of Ref. of the order of 
'one'. This observed feature of Tj/i can clearly be seen from ( fT2l l as 

1=^^^,^{K^)-\ (13) 
s S 

where inverse of kS, in together has an exponent ~ 1.2 for liquid-gas phase transi- 



tion 1151 



For various fluids as discussed in Refs. 11 12112311 one can parametrise 77 /s with the 



reduced temperature, e — \ T /Tq — 1 |, in the domain Tc <T < 2Tc as 

1 



s 



(ae^+b), (14) 



where a and b are different for different fluids. For QGP H, 121 the elliptic flow data 



in RHIC restrict those as 5 ~ 1, a ~ 0.64 - 2/3 and - 0.1 - 1/(4?:), which 
describe the viscous mode despite the various microscopic features in it. Obviously, 
the new fluidity measure J^^i^ in (fT2l) then puts a constraint on the heat flow mode in 
QGP, which in turn restricts the associated transport coefficient K as 



^r^K{ae^ + b)-'^ , (15) 



where the correlation length, ^ is related to various moments of the net baryon number 
fluctuations 

In view of nearly perfect fluidity nature of QGP observed in RHIC we suggest 
through a new exponent-free fluidity that the thermal conductivity, which is a measure 
of transport of energy by the particles, also becomes an important quantity to look at. 
This has not been discussed before in the literature in this context. Interestingly, it has 
been proposed to measure the thermal conductivity in RHIC with its upgradation 1I26I1 . 
which could then test the above estimation of this quantity both in hydrodynamic al and 
nonhydrodynamical regions, and in turn our proposal of new fluidity measure. It will 
also be very interesting, if one can find a possible way to check this universal behaviour 
of fluidity through lattice QCD calculations and in nonperturbative models. The Large 
Hadron Collider (LHC) at CERN will produce a QGP at much higher temperature (at 
least twice of RHIC) and ji close to zero, which would then cool and pass through the 
critical region and it would be interesting to see whether there could still be a Rayleigh 
peak in this limit. 

We note that our analysis is based on the mode-mode coupling theory and the rela- 
tivistic Navier-Stokes theory. A basic assumption of this mode-mode-coupling formu- 
lation is the nondivergence of shear viscosity whereas a proper RG calculation reveals 
that shear viscosity has a very weak divergence governed by the exponent —0.04. On 
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the other hand a Navier-Stokes theory is a kind of mean field equation and may not be 
very appropriate to discuss the dynamics around critical point. The dynamics of the 
fluid in and around the critical point should, however, be described by the evolution 
equation derived from a microscopic theory. Nonetheless, the present understanding 
of the critical dynamics of the produced QCD matter in RHIC experiments is meagre 
and one needs to depend on a kind of mean field description to quaUtatively understand 
the behaviour of the system in and around the critical point because the hydrodynamic 
variables become comparable to the time scale of an order parameter around the critical 
point due to critical slowing down. 
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